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Executive Summary 


This report covers the results of activities of NASA 
Grant NAG3-1106, and extensions: "Analysis of Radiative and 

Phase-Change Phenomena with Application to Space Based 
Thermal Energy Storage". The simplified geometry for thi 
analysis is an infinite, axis symmetric annulus with a 
specified solar flux at the outer radius. The inner radius 
is either adiabatic (modeling Flight Experiment conditions), 
or convective (modeling Solar Dynamic conditions). Liquid 
LiF either contacts the outer wall (modeling ground-based 
testing), or faces a void gap at the outer wall (modeling 
possible Space-based conditions). 

The analysis is presented in three parts representing 
sequential stages of development: Part III, the mitia 

interim report, considers an adiabatic inner wall and 
linearized radiation equations; Part II adds effects of 
convection at the inner wall; and Part I includes the eff 
of the void gap, as well as the previous effects, and 
develops the radiation model further. Although the results 
of Parts II and III are preliminary, and constitute 
background material for Part I, they are nevertheless 
included here for reference, and because they contain 
details not found in Part I, which concurrently with this 
report is submitted as a journal publication. 

The original scope of the grant was to investigate 
analytically the effects of internal radiation upon the 
phase change processes in monocrystal line LiF, with 
extension to effects of polycrystalline structure. However, 
consultations with D. Namkoong during the summer of 1990' 
indicated that the presence of a void gap could have a much 
stronger radiative effect than the change in properties 

represented by a polycrystalline structure. This, indeed, 

turns out to be the case as examination of the results in 
Part I will show. 

The question of the structure of LiF solidifying in 
vacuum under microgravity conditions is moot. Visual 
observation at room temperature of the eutectic, LiF-CaFz , 
solidified in canisters under 1-g conditions, show a 
partially transparent (translucent) optical property for the 
visible spectrum. However, recent experiments and analyses 
of radiation in clear and cloudy ice, indicate a relatively 
small effect of the crystal structure on the overall phase 
change effect. For LiF "cloudy" properties are not 
available and it is felt that using the thin film data o 
Palik and Hunter is the most accurate approach at present. 


The main result from the analysis is the considerable 
differences in melting behavior which can occur between 
ground based 1-g experiments and the microgravity Flight 
Experiments. In the ground based tests, under axial l~g 
conditions, melted PCM will always contact the outer wall 
having the heat flux source, thus providing conductance 
from this source to the phase change front; for this case, 
melting was found to occur primarily from the outer wall > 
with radiative effects causing a small amount of melting 
from the inner wall. In Space based tests and applications 
under microgravity conditions, where a void gap may likely 
form during solidification, the situation is reversed: 
radiation is now the only mode of heat transfer (under the 
axis symmetric conditions of this analysis) and the majority 
of melting takes place from the inner wall. Concurrently 
there is a large temperature excursion in the outer wall 
facing the void gap. However, in both cases, complete 
melting occurs in about the same time, for the adiabatic 
inner wall condition (which is well approximated in the 
Flight Experiments). 

Another major result is the difference between 
adiabatic and non-adiabatic boundary conditions. When there 
is convection at the inner wall, it was found that non- 
melted, partially melted, and fully melted conditions exist 
depending upon the level of the source flux at the outer 
wall. Indeed, the flux level for the planned Flight tests 
would produce no melting at all if the Space Station 
"Freedom" convective conditions were applied. Therefore, 
results from the Flight Experiments can not be extended 
directly to applications with convection, nor can ground 
based test results be directly extended to the microgravity 
environment . 

Overall, the results of this analysis support the 
requirement for interpretive analytical/numerical models m 
conjunction with flight experiments, and it is hoped that 
they provide useful fundamental information and insight. 
However, difficulties associated with void formation and its 
prediction are (perhaps) better addressed with x $dcsi . g n , 
rather than too extensive numerical modeling. For example, 
incorporation of ribbed surfaces on the canister interior 
walls (similar to heat-pipe groves) would provide 
conductance to the PCM, even in the presence of voids; such 
surfaces are currently being investigated under ESA programs 
in Germany. Even more promising is the concept of 
introducing a capillary mesh, or matrix, within the canister 
volume; not only does this provide void control, but it 
enhances the effective PCM conductivity, as well; 
preliminary analyses at UC-San Diego, and experiments at 
local industry, indicate potential for vastly improved 
operation and weight reduction. For future and ° n 9°i n 9 _ 

Solar Dynamic receiver development, it is recommended tha. 
such redesign efforts receive a high priority. 
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ABSTRACT 


A one-dimensional thermal model is developed to evaluate the 
effect of radiation on the phase change of LiF in an annular 
canister under gravitational and microgravitational conditions. 
Specified heat flux at the outer wall of the canister models 
focussed solar flux, or electrically simulated flux; adiabatic 
and convective conditions are considered for the inner wall. 

A two-band radiation model is used for the combined-mode 
heat transfer within the canister, and LiF optical properties 
relate metal surface properties in vacuum to those in LiF. For 
axial gravitational conditions the liquid LiF remains in contact 
with the two bounding walls, whereas a void gap is used at the 
outer wall to model possible micro-gravitational conditions. 

With outer-wall initial conditions at the melting temperature, 
and with the specified flux condition, it is shown that the 
phase-change process is quasi-steady, leading to a simplified, 
but nonlinear system of equations. 

For the adiabatic cases exact integrals are obtained for the 
time required for complete melting of the LiF. Melting was found 
to occur primarily from the outer wall in the 1-g model, whereas 
it occurred primarily from the inner wall in the ti-q model. For 
the convective cases partially melted steady-state conditions, 
and fully melted conditions, are determined to depend on the 
source flux level, with radiation extending the melting times. 

It is concluded from this study that radiation is an important 
effect to include, and that fundamentally different behaviors 
may occur for different gravitational and boundary conditions. 


NOMENCLATURE 


A designates adiabatic wall-2, area [m 2 ] 

B designates Biot-type (convective) wall-2 

C designates liquid contact at wall-1 

c specific heat [J/kg K] 

f radiation function (from Appendix A) 

G designates void gap at wall-1 

g s scaled heat flux source (=g sm /€ 0m ) 
g sm scaled heat flux source (=q s /ffT m 4 ) 

h heat transfer coefficient [W/m 2 K] 

H heat of fusion [J/kg] 

I radiation integral (from Appendix A) 

K scaled linearized radiation conductance 

k conductivity [W/m K] 

N radiation/conduction number (=e om N m) 

Njjj radiation/conduction number (=irT m 3 Ar/k e ^) 

N' radiation/conduction number (=NK a ) 

P radiation exchange function (from Appendix A) 
Q rate of heat transfer [W] 

q s heat flux source [W/m 2 ] 

R thermal resistance [K/W, K 4 /W] 

r radius [m] 

A r PCM overall thickness (=rj-r 2 > 

T absolute temperature [K] 

t time [s] 

t c time constant [s] (=P^,H s ^Ar 2 /k e ^T m ) 

u,v logarithmic functions 

w width (thickness) of wall [m] 


a wall energy storage ratio (=w 1 p w c w T m A () p^H s ^) 

3 Biot number (=h 2 Ar/k e £) 

5j outer scaled liquid thickness (=(rj- r jm)/ Ar ) 

8 2 inner scaled liquid thickness (= (r 2 m~r 2 )/Ar) 

€ emissivity 

<t> flux index 

r radius ratio (=r 2 /rj) 

K conductivity ratio (k es /k e ^) 

A 1 - r = Ar/r x 

n radiation parameter («0.163) 

p density [kg/m 3 ] 

a Stefan-Boltzmann constant (=5.67 x 10” 8 W/m 2 K 4 ) 

r scaled time (=t/t c ) 

r' scaled time for adiabatic wall-2 (=Ng s r=tq s /p^H s ^Ar) 

t" scaled time for conv. wall-2 (=(l-^f)r=tk^(T m -Tf)/p^H s ^Ar 2 ) 
e scaled temperature defect (= (^-1)/Ng s ) 

^ temperature ratio, (=T/Tm) 

Subscripts: 

0 vacuum, initial conditions 

1 wall-1 (outer wall) 

2 wall-2 (inner wall) 

3 surface-3 (at gap facing outer wall) 

a a-band (transparent band) of spectrum 
b b-band (optically thick band) 

e effective 

f fluid boundary 

fully melted 
h convective 




i surface index (= j or 2) 

j surface index (= 1 or 3) 

jc conduction 

l liquid 

n melting condition or temperature 

r radiative 

s solid, source 

st storage 

v wall 


v, page is 
c ; . : quality 


1. INTRODUCTION 


The behavior of phase transition is central to the understanding 
of phase-change thermal energy storage for development of Solar 
Dynamic Space Power [e.g., Labus et al. 1989]. Therefore, ground 
tests have been conducted [Strumpf and Coombs 1990], Space 
Shuttle flight experiments have been planned [Namkoong 1989a, 
1989b] , and numerical analyses have been performed to determine 
two and three dimensional effects of the phase transition process 
[Kerslake and Ibrahim 1990, Wichner et al. 1988]. These analyses 
and experiments utilize an annular canister containing the phase 
change material (PCM) , where a solar heat flux (or electrically 
simulated flux) is impressed on the outer wall of the canister 
(radius r ± in Figure 1) ; the inner wall (radius r?) either is 
convective, or is nearly adiabatic during the heat addition. 

In application to heat receivers the PCM's are high-temperature 
salts, such as lithium-fluoride (LiF) with a melting temperature 
of 1120 K (1556 F) [or the eutectic LiF-CaF2 which melts at 1040 
K (1412 F) ] . At these temperatures radiative transport can be a 
significant part of the overall heat transfer processes within 
the canisters. In the previous two-dimensional model [Kerslake 
and Ibrahim 1990] radiation was not included, and in the 
three-dimensional model [Wichner et al. 1988] the effect of 
radiation is obscured by the complexity of the numerical 
computations; in the flight experiments [Namkoong 1989] 
radiation will occur naturally, but the radiative effect cannot 
be determined directly and must be deduced from canister surface 
temperatures. There is, thus, a need for basic models and 
solutions for the high-temperature enclosure with PCM and 
radiation; to this end a one-dimensional analysis is here 
presented . 

The high-temperature salts exhibit considerable contraction upon 
solidification,' which can lead to void formation at the outer 
wall under microgravity conditions in Space. Hence, in this 
case, radiation is the dominant mode of heat transfer, in 
contrast to axial 1-cf conditions where liquid PCM contacts the 
wall and conduction is dominant. Four cases are considered here 
which model behaviors for adiabatic and convective inner wall 
conditions, and under axial 1-g and n-g conditions. 

Radiation within the canister enclosure follows conventional 
spectral exchange between surfaces, except that the intervening 
medium, LiF, has spectral properties which differ from those of 
vacuum. Therefore, accurate application requires spectral 
integrations; these are presented in Appendix A where metal 
surface properties in vacuum are related to those in the presence 
of the LiF medium. Following Williams [1988], and Song and 
Viskanta [1990], a two-band approximation is made utilizing 
recent measurements of LiF optical properties [Palik and Hunter 
1985] . 

A two-surface thermal model is developed in Section 2 for 
adiabatic and convective boundary conditions at the inner radius. 
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and imposed heat flux at the outer radius. This differs from the 
usual Stefan problem where the temperature at a boundary is 
suddenly changed [e.g., Yao and Prusa 1989, Burmeister 1983, or 
Solomon 1981], resulting in a boundary layer growing in time as 
the error function [e.g., Arpaci and Larsen 1984]. However, as 
shown in Appendix B, with the presently imposed heat flux 
boundary condition, this layer does not develop, with the result 
that PCM temperatures change only slowly in a quasi-steady 
fashion, as was previously observed [Song and Viskanta 1990]. 

For the liquid PCM it is assumed that the conduction limit 
applies, such that natural convection boundary layers do not 
develop significantly and such that the phase-change boundary 
remains axially uniform. This is valid for Rayleigh numbers 
below 1700, for rectangular enclosures [Incropera and DeWitt 
1985, p. 401], but has also been demonstrated experimentally for 
much larger Rayleigh numbers with a heat flux boundary and an 
aspect ratio of 4.5 [Zhang and Bejan 1989]. This is considered 
an excellent approximation under 1-g axial gravitational 
acceleration; under microgravity conditions the magnitude of 
natural convection is suppressed, although some convection 
effects can be present [Bayazitoglu and Lam 1987, Arnold et al. 
1990]. 


2. PROBLEM FORMULATION 

The thermal model is for a solid region surrounded by liquid 
at either or both radii r 3lll and r 2m , as indicated in Figure 1, 
and (in some cases) having a narrow void gap near r±. The gap 
models microgravity conditions where the PCM could completely 
solidify, radially outwards, without making contact with the 
outer wall. Under axial 1-g conditions, there will be contact 
with the wall by liquid filling the annular gap. In these 
annular liquid gaps convection will be driven by the 
gravitational and temperature fields, but without feedback to the 
thermal process for small Rayleigh numbers and long aspect annuli 
[Arpaci and Larsen 1984]. In the following, four cases of 
analysis are considered: A and B referring respectively to 
adiabatic and convective boundary conditions at wall 2, and C and 
G referring respectively to contact and gap conditions at wall 1. 

The heat transfer processes shown in Figure 1 constitute a 
combined radiation and conduction network. Radiation occurs 
between walls 1 and 2 in the transparent spectral band of the PCM 
(the a-band) , and to a lesser extent between surfaces 1 and 3 in 
the optically thick band of the spectrum (the b-band) , as shown 
in Appendix A. In each case considered, the initial condition is 
fully solidified PCM at wall 1 (or surface 3) , and at the melting 
temperature, T™. Under this condition the sensible heat terms 
do not enter the problem, and the quasi-steady network indicated 
in Figure 1 provides the whole solution, as shown in Appendix B. 

With reference to Figure 1, the conduction heat transfer 
from surface i to im (i = j or 2 , j = 1 or 3) is given by: 


Qki = 


( 2 . 1 ) 


T i ~ Tjm 

R ki 

where resistances are 


Rki - , 3-1,3; R k2 - y 1^ .21 

K] 2*k eA J 2*k et 


and the effective liquid conductivity is k e ^. Similarly 

• Ri, = ^( r jm/ r 2 m? 3 

' R *s 2 7r k, ' 3 1/3 


Qks = ^ 


m ~ Tprn 


R ks 


es 


Qh2 * h 2 A 2 (T 2 - T f ) y Q s ti - (pc) w A iWi dT x /dt 
Energy balances yield, respectively, on surface 1 


Qs “ Qra 
on surface 2 


r 0^2 ; C - - contact 
+ ^ Q r b + Qstl • G “ " 9ap 


-< 


Qk2 » A - - adiabatic 


ra \ q^ 2 _ Q^ g • b - - convection 
and on surface 3 

Qrb = Qk3 

The radiative heat rates are related to temperatures as 
Qra — 2Trr^€ om^m^ R a(^l'^2) 

Qrb - 2*^1 *Pb(*i,* 3 > 

where, from Appendix A, 


P a = ^i 4 * 5 " *2*‘° l2(^1^2> 


4.5 


= \b* 4 - 


P b » *1 


+ 3 


(2.2,3) 

(2.4,5) 

(2.6,7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

(2.12) 

(2.13) 

(2.14) 


and u. * 0.163 is a small parameter. 

At the phase boundaries the heat added results in the 
movement of the boundaries through H s ^, the latent heat: 

Qkj “ Qks = “ r jm < j b 1n ' D =1 » 3 

Qk2 + Qks = + 2ir P-6 H s£ r 2m^^ m 


(2.15) 

(2.16) 
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In scaled, nondimensional terms the above relations are 
combined and summarized as follows, 


Wall 1: 

Ng s - NP a + { 

Wall 2; 


(*1 - IJ/v-l , C 
mNP]-, + ctd/'j/d t , 


G 


(2.17) 


±2 - 1 + 

uj ^ 


- *f) = NP a 


_ r Vn{l + A5 2 /r}/A , A 
U 3 ^ ^{(1 - A3j)/r}/KA , B 


Surface im fi = 1 or 3) : 

d«-s - 1 


{! - AA j }- dr - 


j "* I r 0 , A 
Vj " V (1 - 


^ 2 )/ u j » B 


Vj = i*,{l/(l - ASj) )/A 
Surface 2m; 

{r + A8 2 >57 2 - NP a , a 
Surface 3; 


^*5 " 1 

— = * N p b , G 


Here, the radiation/conduction number is 

3 

i 

^e-6 

the Biot number is 


N = lilm^ m 3 (rj - r 2 ) . 

” : 6 Om w m 


(2.18) 

(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 

(2.24) 


a « h 2 (rj ~ r ? ) 

the liquid thickness ratios are 


(2.25) 


8 j 

= yj - r jn , _ £2m-^ 2 

(2.26,27) 

r j “ r 2 2 r j - r 2 

where r 

= t/ , 



P;H c ,(rj - r,) 2 
T m^e^ 

(2.28) 

9s 

_ 2s 2sm 

£ Om 0 ' ^m f 0m 

(2.29) 


4 


a 


( 2 . 30 ) 


P TJ Wl 

r j - r 2 H S^ 

and where r = r 2 /rj , A = 1 - r, and K = k es /k e ^. 



3. RESULTS 


3.1 Case C-A ; 

Under axial 1-g conditions with liquid contacting wall 1 
(nQi»_ case sl . the energy storage in the canister walls is 
negligible and the impJessed . heat flux effectively acts directly 
on the PCM. With an adiabatic wall-2, conservation of energy 
requires all of the energy to go into phase change; this result 
in a closed form solution of the system, (2.17) to (2.32), 
obtained from the sum of (2.20) and (2.23): 


- AS ! 2 /2 + rs 2 + AS 2 2 /2 * Ng s r = r' 

In particular, for complete melting when r 2m = r im ( 5 l + 
the scaled time is obtained from (3.1) as 

l»9s T fm “ r' tm - (1 + r)/2 

or, in dimensional variables, the fully melted time is 

2 - r, 2 )o,lU, 


t fm 


„ *(ry 


(3.1) 

= 1), 

(3.2) 

(3.3) 


27rr 1 q s 

Thus, for case C-A the time to melt the PCM is the Phase-change 
heat content, divided by the input heat transfer rate. 

Because of the liquid contact with the canister walls, the 
wall temperatures do not depart strongly from T ro ; therefore, for 
the "C" cases, it is adequate to use the linearized radiation 
exchange, P a (*i,* 2 ) - *a<*l ” *2> » . where ' fc 5 e constant is 
K a * 4.5 f 0 c(l), and f 0c is from Figure A-2 of Appendix A. 

For the parallel plate limit (A > 0, r “7 _> i) the 

adiabatic, linearized system has the exact solutions: 


r' 

2 + N'r'l 

2 

1 + N'r'J 

r ' 

2 + N'r'l 

4 

1 + N'r'J 

r' 

2 + N't'J 

2 

1 + N'r'J 


= t* f N 1 r 1 
2 [l + N'r ' 


(3.4,5) 


(3.6) 


( 3 . 7 ) 


where N' = NK a , r 


= Ng s r and 6 - (^ - 1)/Ng s . 


These solutions clearly show the double surface phase-change 
phenomenon (as also indicated in Figure 2-1, A) • *t is seen in ^ 
equations (3.6,7) that Si grows linearly with r for small r 
(i.e., initially), whereas S 2 grows quadratically ; the initially 
linear phase boundary growth has been observed previously in the 
absence of radiation [Evans et al. 1950]. 
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Wall temperatures are determined in this analysis, rather 
than specified. The gradual increase in the wall temperatures in 
(3.4,5) shows the absence of thermal boundary layers, as proved 
in Appendix B. From (3.4) the maximum wall temperature increase 
at complete melting when r' = 1 is /'i max ” 3. = N 9s ; i- n PkY s i ca ^ 
terms T-, max - T m - q s (ri-r 2 )/W, or about 25 K above 1120 K, 
which validates linearization of the radiation term. 


In the case of an annulus with finite radii (r 2 /rn = r * 1) , 
equations (2.17) to (2.22) were solved by numerical integration, 
with results as shown in Figure 2, for r = 0.5. The initially 
linear and quadratic growth of the two phase fronts is evident; 
complete melting occurred at = 0.75, as predicted by (3.2) . 


3.2 Case C-B ; 

With the convection heat sink at the inner wall, the 
possibility exists that q s may not be sufficient to cause any 
melting at all. This is seen from the combination of (2.17) to 
(2.20) which yields 

(1 - ASi)^ 1 = Ng s - 3 T(i/' 2 “ *f) (3.8) 


where the initial slope must be positive for to grow; thus, 

Ng s must be greater than the minimum 

Ng s0 « s r<* 20 - *t) = sr n , N 1 I/up 0 ; - »f> < 3 - 9 > 

where from (2.19) u 0 = }/KK , and * 2 p = *2(°) is obtained 

from (2.18). In dimensional variables with h rl = K a' the 

limiting minimum heat flux for phase change to occur is 


QsO “ 


r i 

h 2 r 2 


Tm - Tf 

+ ^ 

h rl r l + k s /^{r 1 /r 2 ) 


(3.10) 


which is the initial overall temperature difference over the 
total thermal resistance. Similarly, there is a critical heat 
flux for which the PCM just becomes 100% melted and at steady 
state with zero slope at 3^ = 3; this is given by (3.8) as: 


Ng sl - 8T(1 - * f ) , or, q sl = h 2 (T m - T f )r 2 /r! (3.11,12) 

Partially melted, and unmelted, conditions have also been 
observed experimentally [Strumpf and Coombs 1990]. Finally, 
there is a limiting flux above which melting occurs from both 
sides. 
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•m 




Ng s2 - Ng sl 


1 + 


(3.13) 



however, for small N', this is a high flux which is usually not 
encountered. 

To summarize, g- < g s0 results in no phase change, g s0 < g s 
< 9sl results in a steady state with only partial melting, g sl < 
9s < 9 s 2 results in 100% melting in finite time from wall-1, and 
9s — 9s2 results in 100% melting in finite time from both walls. 

The numerical integration of the system equations for the 
partially melted case is shown in Figure 3, where the scaled heat 
flux was taken as 


9 S ■ *9 S 1 + (1 - *)9 S 0 ' 0 £ * < 1 (3.14) 

For the parallel plate limit the partially melted steady-state 
condition can be obtained analytically as 

8 m N* + k + g 

!,ss * N' + K + 03 + (1-0) (K-IJN' ( 3 * 15 ) 

For ?s > 9sl» the pase-change front grows rapidly to complete 
melting, as shown in Figure 4, where the applied flux was taken 
as g s = 0g sl < g s2 . 


3.3 Case G-A: 

Under microgravity conditions where a void gap may form 
between the PCM and the outer wall ("G M -cases) . there is onlv 
radiative transport between the outer wall and the PCM. This 
will cause considerable temperature increase in wall-1, and 
some energy storage in the wall. 

There are no further simplifications possible to the system 
(2.17) to (2.23). Nevertheless, for an adiabatic wall-2, it 
possesses the exact integral: 

a (*l - 1) + «3 - AS 3 2 /2 + rs 2 + AS 2 2 /2 = Ng s T = r' (3.16) 

which, for the fully melted condition, S 2 + S 3 = l, reduces to 

T ’fm = (1 + r)/2 + a(* 1/fm - l) (3.17) 

By comparison with (3.1) and (3.2), these results represent the 
sum of the energies stored in the PCM and in wall-1. Exact 
calculation of the melting time from (3.17) is not immediate 
because the original system must be solved to determine wall 
temperature, at any time; however, the ratio of wall 
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Ca E a «/ tance to Phase-change capacitance is a small quantity, 

* “ so that melting times are only slightly longer than 

for liquid contacting wall-1. 

A numerical integration of system (2.17) to (2.23) for 
r 0.5, g sm< — 0.1, a— 0.1, and for radiation functions from 

sllowl1 in Figure 5. In contrast to case C-A, most 
of the melting takes place from the inner wall, and only a slight 
near the outer wall; otherwise the melting is completed in 
about the same time as previously because of the small value of 
a. 


The corresponding wall temperature variations are shown in 
Figure 6 , where there is a dramatic sudden increase in the outer 
temperature, in contrast to the liquid contact case where 
this temperature increased only slightly and gradually. 
Asymptotic analysis shows this initial temperature "jump” to be 
approximated by J F 

*l,outer(°) ~ C 1 + 9sm/ e Omfogf 1 ) j 1/ ' 4 ' 5 (3.18) 

The effect of the emissivity on is evident in both (3.18) and 
Figure 6, whereas there is a smaller effect on the interface 
temperature, * 3 , increased only slightly from 1. Clearly from 
(3.18), an increase in the applied flux, g sm , would cause a 
further increase in the outer wall temperature. 


3.4 Case G-B; 

With the gap at the outer wall, and convection at the inner 
wall, the flux limit definitions for g s0 , (3.9), and g c1 (3.11) 
apply as previously; however, * 2 must now be calculate! from the* 
nonlinear equations. Because v. is a small quantity it is found 
for the steady -state limit that *20 = 1 “ °(^N), such that g s0 = 
”sl(-*- ” ®(*0); that is, for any phase change to occur at all, g_ 
must be very near the limit, g sl . Above this limit, g s2 = g s1 {! 

+ , ?(^)>» so that melting will occur from wall-2 for g-iust 1 
slightly above g sl . * s J 

An example of this behavior is shown in Figure 7, where g_ = 
9si: ,£?r & — 1 there is only slight melting from surface 3, and 
no melting at surface 2; for * = 1.5 the melting rates are about 
the same from both surfaces, but a long time is required for 
complete melting; for <t> > 2 the melting rate is greater from 
surface 2, the inner wall, and complete melting was attained for 
times, r , shown in the figure. Comparison with the 
adiabatic cases in Figure 4 shows that a longer melting time is 
needed with convective heat transfer, as expected. 

. Corresponding temperatures are shown in Figure 8 ; after the 
initial Dump, these remain essentially constant during the 
melting process. With the initial condition of 83 ( 0 ) = 0 

/t 3 io( - H' greater than 1, approximately as given by 

(3.18), and * 2 t0) is slightly less than 1, as shown. 
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CONCLUSION 


Four analytical cases have been considered for the melting 
of Lithium-Fluoride in an annulus with impressed heat flux at one 
boundary, including the effect of internal radiation heat 
transfer. It was found that this process is quasi-steady when 
the solid LiF near the outer wall is initially at the melting 
temperature. Radiation was found to be an important effect, 
especially in the presence of void gaps near the outer wall. 

For the adiabatic inner wall condition, the time for 
complete melting is a fixed quantity which depends only slightly 
on voids. However, the location of the phase boundaries is 
strongly influenced by the void gap, with melting occurring 
primarily form the outer wall when there is liquid contact, and 
primarily from the inner wall when there is a void gap at wall 1. 
Wall temperatures remained close to the melting temperature with 
liquid PCM contacting wall 1? but, the wall-1 temperature has a 
large and sudden increase above the melting temperature when a 
void gap is present at wall-1. 

For the convective inner wall conditions there may, or may 
not, be complete melting, depending on the level of impressed 
heat flux, relative to the fluid temperature and other problem 
parameters. For liquid PCM contacting wall 1, melting occurred 
only from wall 1 for moderate heat fluxes, and steady-state 
partially melted conditions were determined; for a void cjap at 
wall -1, conditions for partial melting almost did not exist, and 
a substantially high heat flux is required in order to cause any 
melting at all. 

In summary, the results show fundamentally different 
behaviors depending on the boundary conditions employed. Thus, 
experimental results from the flight experiments, where the inner 
wall is essentially adiabatic, can not be used directly for 
convection applications, without interpretive 
analytical/numerical modeling; nor can ground-based results, 
under 1-g acceleration, be used directly in microgravity 
applications. 
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APPENDIX A 


RADIATION MODEL 


Experimental data for Lithium Fluoride indicate a near 
perfect transparency for wavelengths below 5.5 Aim, and optically 
thick properties for wavelengths above about 7 *m [Palik and 
Hunter 1985]. Therefore, a two-band radiation model is 
considered where the LiF is transparent for wavelengths below 
5.5 irn (the a-band) , and optically thick for wavelengths above 
this value (the b-band) . (This model obviously neglects 
processes in a narrow band where the LiF is neither transparent 
nor optically thick. ) 


In the a-band there is radiation exchange between the two 
bounding metal walls, and spectral variations of properties of 
the walls, and of the LiF with index of refraction n, are 
included. In the vicinity of T- = 1120 K most of the emitted 
radiation is in this transparent band. In the b-band, there is 
a minor radiative addition to the LiF conductivity and, in the 
prescence of a void, minor radiation between wall-1 and the LiF 
surface. 


Wall Emissivitv: 


Spectral emission in to a medium with refractive index n is 
given by the Planck formula and the spectral emissivity [e.g. 
Sparrow and Cess 1978]: 


e v ~ ^v e bv n ^v e bc — e r e bC 

where C = c 2 /nXT = v/v 0 , v 0 = c 0 T/c 2 , c 0 = 3xl0 8 m/s, 
c 2 = 14,388 Aim K, and where 

15 <jT 4 ? 3 
7r^ v o e^ - 1 


(A.l) 


(A. 2) 


Here € v represents emission into the medium relative to blackbody 
emission into the medium, and e* represents emission into the 
medium relative to blackbody emission into vacuum : the latter can 
be greater than 1 whenever n is greater than 1. 


For a number of metal _ surf aces at the high temperatures 
under consideration, experiments have shown that the normal 
spectral emissivity varies with wavelength to the -1/2 power, 
even into the visible part of the spectrum [Seban 1965, 
Touloukian 1970] ; therefore, consistent with electromagnetic 
theory, we take 




e X 


n € xc 



M 


nx 


n e Xc 



(A. 3) 


where e Xc is the best-fit emissivity for n = 1 at wavelength X c , 
and where C c = c 2 /X c T. For emission into vacuum this model 


1 


yields the total normal emissivity 


6 0t = 


■ft? 


15 r" t 3 - 5 ac _ 1_ t,„ 

** J 0 eC • 1 c * NfT7 


(A. 4) 


where c 4 = ir 4 / (15x12 . 27) * 0.529. For example, for Nickel at 
1390 K with the experimental value e^c = °* 25 a ^ *c = 1 C Se k an 
1965], the model yields € nt = °* 147 ' f° r i \c ~ °*°® at Xg = 9 
it yields 6 ot = 0.141, which indicates good agreement with 
experimental measurements for this material at this temperature. 


Now combining (A.l), (A. 2) and (A. 4), the effective spectral 

emissivity for emission into medium n is given by 


= c 4 n 3 e ot^ f (A. 5) 

Thus, if the total normal (hemispherical) emissivity function 
£ 0 t(T/T m ) for emission into vacuum is known in the vicinity of 
Tjj 7 then (A. 5) yields an approximate spectral noraal 
(hemispherical) emissivity for emission into medium n. For 
example, with € = 0.147 from the above calculation, the 
emissivity at X = 1 vrn for emission into LiF (n * 1.4), is 0.58. 


A-Band Radiation Exchange; 

Application of radios ity microbalance to the phases and 
surfaces indicated in Figure 1 yields the spectral flux at wall 
1 [e.g. Siegel and Howell 1981]: 


sri = 


e ^i ~ *hl2 

+ iZIj_2 _ Te 

£ C1 € C2 


(A. 6) 


where the effective interface transmittance is 

T e = (l - 2p 3 ) / ( 1 - p 3 ) » r 3 * 0.98. Thus, the a-band total 

radiative flux is 


<3al = J! Svl^v = v 0 J" q cl dC 


Combining (A. 2) with (A. 5) to (A. 7) results in 


T = ^1 4 * 5 Ii - $2 


6 0m a T m 


4 - 5 I 2 


(A. 7) 


(A. 8) 


where 


VW.iK 
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I-i = 


12.27 


fe^ ~ 1) dc 

„nT77 


n l 3 € Ot(^l) F 12 n 2 3 € 0t(^2) 


n \T7r 


■ c 4 € 0m T 3^ *i f ' 


(A. 9) 


an d C a i = Fam/^i' ^am - c 2 v a/ c 0 T m' “ T i/ T m' 6 0m “ ^Ot^ 1 )* 

In (A. 9), n^ = nf^C) and n 2 *= nf^f) for LiF contact with 
both walls, whereas n 1 = 1 in the prescence of a void gap at 
wall-1; for the region considered, the LiF spectral refraction 
index data of Palik and Hunter [1985] may be represented as a 
function of wave number, as follows: 


n(iK) = 1.38 + 7.79X10 -4 *r - 2.3 e" 1 * 87 ** 


(A. 10) 


A closed-form integration of (A. 9) is not possible; however, 
it may be evaluated in the vicinity of T m by use of Taylor series 
expansions about this temperature. Since wall temperatures, on 
an absolute scale, do not depart strongly from T m , this is a 
valid procedure, which yields for LiF contact with wall 1 and 2, 

I i,c = ^Oc(^i) 

+ Oi-i + ^r i ]{(7^ n - j]fic<*i--n(n) + f 2 c (n>} 


+ € Om 

and for a void gap at wall 1, 


Ii,g “ t 0 g(*O + 

+ ' 5J f ig^i' n (n) + f 2 g(»i>} 


+ £ 0m r 3 (*i-l)f3g( + i) 

j 

where for 7 = c or g, j 


f 07^i) - 12 ] 27 j 


^am/^i d 07 




(A. 11) 


(A. 12) 


(A. 13) 


(A. 14) 


f 27^i) 12.27 I 
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f , = c 4 / 2 _ [14 

37( * lJ 12.27 Jfam/^i 

Here the denominators are 


NfPqfndf 


'07 


Doc - 1 - =4‘on^T 


D 0 g - 1 + 


F 12 i 3 (C) ‘ ° 4 £ 0 ^ 3 'Tr 


and the numerator function is 


g(f) 


? 3 • 5 


(A. 17) 
(A. 18) 


(A. 19) 


which is less than 1% at the computational upper limit of f = 14; 
here the apostrophe ( * ) denotes the derivative with respect to 
the function argument. A sample evaluation of f 07 is shown in 
Fig. A-2 . 


B-Band Radiation Effect; 

In the optically thick region of the spectrum, there is 
radiation exchange across the wall-1 void gap according to 


<3bl = 



e hf1 ~ e hC3 - dr 

— ~ + — 1 

e Cl £ C3 


(A. 20) 


where is obtained as before, but with n ** 1, and €^ 2 I s 

obtained from electromagnetic theory as 


_ 4n 

€ t 2 “ (n + 1) 2 + k 2 


(A. 21) 


With the data of Palik and Hunter [1985], evaluation of (A. 20) 
and (A. 21) resulted in 

<3bl = 0m ffT m 4 (^l 4 - *3 4 ) (A- 22 ) 

where n is a slightly decreasing function of increasing 
temperature, with an average value of n - 0.163. Thus, radiation 
exchange across the gap is quite small. 


Energy absorption in the medium is modeled as Rosseland 
conductivity [e.g., Edwards 1981, p. 296; Ozisik 1973, p. 318]; 


(A. 23) 


*R 



_i_ , 9e hi 
3a x dT 


dx 


5 



C 4 eC 

a x (e^ - 1)2 


d? 


where a x = 4irk/X . Evaluation of this function showed but a weak 
temperature dependence, and an approximate value of 0.005 W/cm K, 
which is an order of magnitude less than the medium thermal 
conductivity. Thus, there is only minor b-band radiative effect 
on the phase-change heat transfer problem. 



Fig. A-l. Radiation Function Evaluation. 
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interval , 0 < t < t 1( the wall temperature 
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singular perturbation problem in time, which may be solved by 
asymptotic expansions or multivariable methods [Nay f eh 1981, 
Smith 1985]; in combination with the integral approximation 
[Arpachi and Larsen 1984] this yields to lowest order the 
composite expansion 


T w - T „1 + + fjf [l - e' 3a(t ' t l ,/L ° 2 ] 


(B. 11) 


It is seen in (B.12) that there is a discontinuity in wall 
temperature as represented by the "inner-time" exponential. But, 
again, for Lq — > 0, this singularity is removed and only the 
outer, "quasi-steady" solution remains; 

q2 

T w, outer = T m + t (B.12) 


It is noteworthy _ that the linear outer-time solution satisfies 
the initial condition Ty(0) = T m , in contrast to the conventional 
temperature Stefan problem where the initial wall temperature at 
tn is not T m . 


Therefore, for the solid PCM initially contacting the 
canister wall, or for small initial liquid layers, there is no 
significant boundary layer effect, or effect of the liquid 
specific heat, and the outer-time solution provides the whole 
solution to the problem. This conclusion is also reached when 
the initial temperature profile is different from constant at T w . 
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PART I I 


RADIATION AND PHASE-CHANGE 


OF MONOCRY STALL I NE LITHIUM-FLUORIDE 
IN AN ANNULAR ENCLOSURE 
WITH SPECIFIED HEAT FLUX 


by 

Kurt 0. Lund 
135 Sixth Street 
Del Mar, CA 92014 


Abstract 


A two-band radiation model is utilized for radiation exchange of 
LiF contained within an annular canister, and its effect on phase 
change is evaluated. Linearization and the quasi steady 
approximation leads to simplified one-dimensional heat transfer 
models, having a specified heat flux at one boundary and either 
an adiabatic or convective condition at the other, and to several 
exact solutions. 

For the adiabatic case, radiation caused phase-change to occur 
from both boundaries, and eliminated the influence of solid 
conductivity; for the convective case, partially melted 
steady-state conditions, and fully melted conditions were 
determined to depend on the heat flux level, with radiation 
extending the melting times. 

Radiative exchange between the two vails was found to be about 
twice as high for the LiF medium as for vacuum, but the overall 
effect of radiation was limited to about ten percent of the total 
heat transfer rate for practical heat flux levels and annulus 
dimensions. Further analysis is required to access the effect of 
voids within the enclosure, where radiation may have a stronger 
effect. 


QB 

or 


I'liWAL 


NOMENCLATURE 


a bsorption coefficienty l/(cs) 

a,b fitted constants 

C modified blackbody radiosity (W/cm 2 ) 

c radiation constant 

e emmissive power (W/cm ) 

f emissivity fractional function 

g scaled, nondimensional heat flux) 

H latent heat of fusion (J/g) 

h heat transfer coefficient (W/cm K) 

k extinction coefficient, thermal conductivity ( W/cmK ) 

N radiation number 

n refractive index 

Q heat transfer rate per unit length (W/cm) 

q heat flux (W/cm 2 ) 

q+ radiosity ( W/cm ) 

R radiation or thermal network resistance (cm 2 ), ( cm-K/W ) 

r electrical resistivity (Ohm-cm), radius (cm) 

T absolute temperature (K) 

t time minutess) 

u,v,w logarithmic conduction terms 


Greek 

a 

£■ 

5 

t, 

r 

e 

K 

■x 

A 

P 

a 

-r 

t 


absorptivity 
convection Biot number 
scaled liquid thickness 
emissivity 

fraction of limiting heat fluxes 
radius ratio 

nondimensional surface temperature 

ratio of ef f . solid to eff. liquid conductivity 

wave length 

1 - r 

density (g/cm 3 ), reflectivity 
Stef an-Boltzmann constant ( W/cm 2 K 4 ) 
nondimensional time 
integration variable 


Subscripts 

0 pertaining to Tq 

a pertaining to transparent upper limit, 

b blackbody 

e effective 


1 


f fluid of convection boundary 

h convection boundary 

i surface i 

k conduction 

£• liquid 

m melting or phase-change temperature 

n normal direction 

r radiation, reference 

s solid, source 

v vacuum 

X spectral quantity 

1,2 surfaces 1 and 2 

0,1,2 flux limits 



1 . INTRODUCTION 


The behavior of phase transition is central to the 
understanding of phase-change thermal energy storage for 

1989l° Pm TherpL SOlar Dyn * mic Space Power C©-8.. Labus et al . 
19 T 9 ^ ' Therefore, ground tests have been conducted [Strumof 

Dl«nS°H m ^M 1 ? 9 ° ] ’ SPaC ® Shuttle flight experiments have been 
bi^n i amk °°r 9 i 989a > 1989b], and numerical analyses have 
been performed to determine two and three dimensional 
effects of the phase transition process [Kerslake and 
Ibrahim 1990 Wichner et al . 1988]. These analyses and 
experiments utilize an annular canister containing the phase 
change material (PCM), where a solar heat flux (or 

of e ?hI 1 ^ al - y * Sim J ,la 5* d flU V is * mpr essed on the outer wall 
,^ h f c ^ lster (radius r i ) , and convection at the inner 

i 2) ' th ? Flight Experiments the convection 
is Simulated by conduction in a rod and subsequent 
radiation to space [Lund 1991] ; the effect is a near 
adiabatic condition at n during application of the heat 

■ 1UX at r . 

In application to heat receivers the PCM’s are high- 
temperature salts, such as lithium-fluoride (LiF) with a 
rSV"? t^Perature of 1120 K (1556 F), or the eutectic 
LiF-CaF 2 which melts at 1040 K (1412 F). At these 
temperatures radiative transport may be a significant part 
of the overall heat transfer processes within the canisters 

l^OWpHfpMo 8 two_dl mensional model [Kerslake and Ibrahim 
1990] radiation was not included, and in the three- 
dimensional model [Wichner et al . 1988] the effect of 
radiation is obscured by the complexity of the numerical 
ln the fll9ht experiments [Namkoong 1989] 
k P w 1 ^ 1 ° PCU L *]? turall y- but the radiative effect 
tlmnpr*r e deter .^ n * d directly as only canister surface 
temperatures will be measured. Combined convection and 

1 i 1 was J f nv ? sti 9fted for a square enclosure 

[ et al. 1989] , radiation with scattering in an annulus 

was considered [Tsai and Ozisik 1990], and plane one- 
dimensional combined radiation, conduction and phase change 
was investigated for low-temperature ablation [Yuen and 
Khatami 1990] and ice removal [Song and Viskanta 1990] 

9 nee ? f ° r b3Sic models and solutions for 
the high temperature enclosure with PCM and radiation; to 

this end a one-dimensional analysis is here presented. The 
present analysis is limited to monocrystalline salts in 
contact with the canister walls; the presence of voids 
and/or polycrystalline structure will be treated in a 
separate publication. 

and ac^? U hLJ h ?, CaniSterS u inVariably are finite in length, 

-S 1 h fluxes may be circumferentially and axially 
n uniform, there are a number of reasons why a one- 
dimensional approximation is useful: 


a) closed-form or simplified solutions are possible 
which illustrate the basic phenomenological 
interactions ; 

b) solutions provide bases against which complex 
numerical models may be compared for identical 
boundary conditions; 

c) solutions provide ready estimation of experimental 
behaviors, such as surface temperature-time 
variations ; 

d) the phase— change process is predominantly one- 
dimensional in the radial direction for canisters 
with larger length— to— thickness aspect ratios. 

This was the case even at a lower ratio [Kerslake 
and Ibrahim 1990] . 

Radiation within the canister enclosure follows conventional 
exchange between diffuse surfaces, except that the 
intervening medium, LiF , has spectral properties which 
differ from those of vacuum. Therefore, accurate 
application requires spectral integration. Following 
Williams [1988] , and Song and Viskanta [1990] , a two-band 
approximation is made utilizing recent measurements of LiF 
optical properties [Palik and Hunter 1985]. It is shown in 
Section 2 that at the melting temperature, Tm = 1120 K, the 
majority of radiation emissive power occurs in the 
transparent part of spectrum for LiF. For the smaller part 
at longer wavelengths, a "thick-gas" approximation is used 
which adds a term to the apparent thermal conductivity. 

The emissive properties of the metal canister surfaces is 
assumed to follow the Hagen-Rubens relation. With the 
spectral properties of LiF this leads to a spectrally 
integrated emissivity fractional function which is the ratio 
of the average emissivity of the surface contacting LiF to 
the emissivity of the surface in vacuum at the same 
temperature . 

A two-surface thermal model is developed in Section 3 for 
adiabatic and convective boundary conditions at the inner 
radius, and imposed heat flux at the outer radius. This 
differs from the usual Stefan problem where the temperature 
at a boundary is suddenly changed [e.g., Yao and Prusa 1989, 
Burmeister _ 1983 , or Solomon 1981], resulting in a boundary 
layer growing in time as the error function [e.g., Arpaci 
and Larsen 1984]. However, with the presently imposed heat 
flux boundary condition, this layer does not develop, with 
the result that surface temperatures change only slowly in a 
quasi-steady fashion, as was previously observed [Song and 
Viskanta 1990] ; this phenomenon is utilized presently as the 
quasi-steady approximation. 

Besides the quasi-steady approximation in the liquid, it is 
further assumed that the conduction limit applies, such that 
natural convection boundary layers do not develop 


“-sttaKSSSSSS® 

ccn Sf ,ered an .xd 5ST? J!f iS 

SSSsWSr^ 

1987, Arnold et al . 1990]. P res «nt [Bayazi toglu and Lam 

t^knts an 

independent of radiative and conductive effects with? i%K 

»a"n J^d tL d vot ndS on l y , on the h« at flux at the r 
the volume end heat of fusion of the PCM The 

energy an'd oaus^^e^^^^?^^^^ H ^-^ibute 
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2. TWO-BAND RADIATION MODEL 


The effect of radiation within the enclosure is determined using 
the optical properties of LiF: refractive index, n, and 
absorptive index, or extinction coefficient, k. A summary of 
recent measurements is shown in figures 2-1 and 2-2 for the 
spectral region of interest [Palik and Hunter 1985]; also shown 
is the Planck spectral emissive power for the melting 
temperature, 1120 K, which shows that emissions whithin the 
canister lie primarily between 1 /.on and 6 ion. It is seen that 
there is a transparent region of the spectrum near 1 Aim (where 
k — > 0); the refraction index appears relatively constant in 
this region, but decreases somewhat with increasing wave length. 
In addition to the optical properties of the LiF, the canister 
surface properties are required, as will be shown. 

2.1 Emissivity and Absorptivity 

The assumption is made in this analysis that the metal surfaces 
bounding the LiF are diffuse. For spectral emissivity, the 

emissive power of the surface is given by [Sparrow and Cess 1978] 
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( 2.1 ) 


where L = c 2 /nVr, L a = c 2 /n \ T = c 3 /T = 1995(K)/T, n* = dn/d*. 
and where c 2 is the radiation constant, c 2 = 14,388 AtfnK . Here, 
the index of refraction for LiF is the spectral quantity shown in 
figure 2-1, which is approximated as 


n = 1.39 - 0.00271-X 2 = a - b* 2 


( 2 . 2 ) 


and the upper limit of integration is taken as > = 5.5 ton. 

For the metal surfaces bounding the LiF the Hagen-Rubens relation 
is assumed for the spectral normal emissivity [e.g.. Sparrow and 
Cess 1978] 

£ Xn = 36 . 5 '4 r / a ' (2.3) 

where r is the electrical resistivity in Ohm-cm, and \ is the 
wavelength in Aim. The result in vacuum of spectrally integrating 


(2.3) is the total normal emissivity 


t, 


n , v 



(2.4) 


where T s is the surface temperature in degrees Kelvin. The form 
(2.3) applies well to resistive alloys, such as stainless steel 
[Edwards and Bayard de Volo 1975, as quoted in Edwards 1981], and 
is assumed here to represent, approximately, other canister 
materials such as the Haynes alloy [Strumpf and Coombs 1990] 


Let r be proportional to the absolute temperature, then r = 
r 0 T s /T 0’ w ^ ere r 0 is the resistivity at temperature Tq , and where 
from (2.4) the total normal vacuum emissivity at this temperature 
is 6 0,v Substitution of these variables into (2.3) results in 
the normal (or hemispherical) emissivity ratio 



36.5 

0.576T 0 


\T 


v* 



(2.5) 


where c^ is the nondimensional number, c 4 = 36 .5/0 .576\| c 2 ' = 
0.529. Now, substitution of (2.5) into (2.1) results in 


e a = * T „(V T o] e 0,v f ca <V 


( 2 . 6 ) 


where the emissivity fractional function is defined by 






n + Xn ’ 
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(2.7) 


This function is shown in figure 2-3 for * a = 5.5 m, where it is 
seen that the effective total emissivity is roughly twice that of 
the value in vacuum: 
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(2.8) 


A completely similar procedure yields the effective absorptivity 


«a< T s’V = £ 0,v 
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f Ca<T e ) 


(2.9) 
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where T 0 is the environment (other surface) temperature. 


For subsequent calculations, the emissivity function is 
approximated as 

f £a - 2-°mn m ) 0 - 278 = (2.10) 

Although the preceding model is approximate, it nevertheless 
incorporates the essential features of emission and absorption of 
the surfaces bounding the LiF in its transparent region. In any 
event, results cannot be more accurate than is known. In 

the subsequent calculations, T Q is taken as T m ’ such that = 

& m,v ancJ t ^ )e a bove temperature dependence need only be accurate 
near the melting temperature, T . 

2.2 Thick Gas Approximation 

It would appear from figure 2-2 that k is sufficiently small for 
the wavelengths of interest that the medium can be considered 
perfectly transparent, there is, however, a region at the longer 
wavelengths where absorption occurs in the medium, as can be 
shown by evaluating the absorption coefficient [Siegel and Howell 
1981, p. 427] from the data of Palik and Hunter [1985]: 


4-JTk 



as shown in figure 2-4. These values may also be read as the 
optical thickness for a 1 cm layer of LiF (or as half the optical 
thickness for a 2 cm layer). The results are in general 
agreement with previous measurements [Amr . Inst. Phys . Handbook 
1957]; thus, the assumption of an optically thick gas above \ = 6 
or 7 Atm is supported by the more recent measurements. The net 
effect is an additive term to the LiF thermal conductivity 
[Williams 1988] : 

16crn 2 T 3 T( K ) - 320 

k _ % 0.02 [W/cmK] (2.11a) 

r 3a, 800 

A 


Since, as seen in figure 2-1, there is but little radiant energy 
at the higher wavelengths, this two-band radiation model is 
considered an adequate approximation, especially in view of 
uncertainties associated with the canister surface properties. 
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2.3 Radiant Energy Exchange 


For each surface, i, in the enclosure the net radiative heat 
transfer rate can be expressed as follows. 


C - q 
i i 

Q = (2.12) 

X R 

i 


where q^ + is the radiosity, where C- = €- crT^/cx^ is the modified 
blackbody radiosity, and where the surface resistance is 


1 - a 

i 

R = (2.13) 

i Cf. A 

i i 


These relations may be combined with formulas for radiosity 
between surfaces [e.g., Edwards 1981, p. 117] so that 
conventional radiation network solutions may be used with the 
modified blackbody surface radiation, and with the surface 
resistance given by (2.13). A closed-form analytical solution 
is, of course, not possible as depends on the "other-surface 
temperatures, as shown in (2.9); however, an adequate iterative 
solution can be obtained with a weighted average for T . For a 
two-surface enclosure, (2.12) and (2.13) yield 


C - C 
1 2 

Q = - Q = (2.14) 

1 2 1 - a 1 - a 

11 2 

+ + 

A Cf. A F A a 

11 1 12 22 


or, with F 21 * 1 , V* - T/T m> and ( 2 - 8 ) to (2.10), the radiation 
exchange is 


°1 - V JT p< w 
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(2.15) 


where 
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(2.16) 
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A linearization of (2. l\> about T 
P( W = K rl (^ - iP 2 ) 


~ T m ~ 1 ) results 
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(2.17) 


where 
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5 + 2s 


rl 1 r 

— [i + 'V« 2 J - i 

m ,v £m 


(2.18) 


Since, as will be shown, wall temperatures differ only 
slightly from T m , the linearization is valid and within the 
accuracy of other approximations made 



3. HEAT TRANSFER MODEL 


The PCM contained by the annular canister in the Flight 
Experiment [Namkoong 1989] , as well as in applications with 
annular geometry [Strumpf and Coombs 1990] , is to cycle between 
the charging (heat addition) and discharging (heat removal) 
modes. In the first mode heat is added to (impressed upon) the 
large-diameter, outer wall to melt the PCM while some heat 
transfer may or may not occur at the small-diameter, inner wall; 
in the second mode, heat is transferred from the inner wall to 
solidify the PCM while the outside wall is, more or less, 
insulated . 

We consided the basic charging mode where at time zero the 
external heat source is Q s ( W/m ) and the PCM is 100% solid and at 
temperature, T m , at the outer wall. Furthermore, to obtain a 
basic (or fundamental) solution, an infinitely long cylinder is 
assumed together with axis-symmetric heat input. Additionally, 
axial gravitational acceleration is assumed such that any void is 
removed to the end region, and such that liquid remains 
continuously in contact with the two walls and the phase-change 
boundaries . 

The above assumptions result in a radially one-dimensional 
time-dependent system, as illustrated in figure 3-1. Because of the 
transparent part of the LiF spectrum there is radiation from the 
outer wall (surface 1 ) to the inner wall (surface 2), represented 
by resistance, R r . Two boundary conditions are considered on the 
inner wall: a) adiabatic condition, which is representative of 
the Flight Experiment [Namkoong 1989], and b) convective 
condition , which represents the Space Station [Strumpf and Coombs 
1990] (although this application uses the eutectic LiF-CaF 2 as 
PCM, the preceding optical properties of LiF are used in sample 
calculations, as seen in Table 3-1). 

In case a the inner wall acts as a refractory or reradiation 
surface, such that the energy source, Q s , causes phase change at 
both melt-surfaces, lm and 2m, as shown in figure 3-l(a). The 
result is that the solid PCM remains at the initial uniform 
temperature, T m , until complete melting has occured; thus, solid 
properties do not enter in to this heat transfer problem 
(i.e., R ks = co). 

In case b the inner wall act as a heat sink whenever T 2 > T.^; 
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therefore. T-, will initially be less than 1,(0) = T 
igure 3-l(b). Thus, no initial melting ocLrs at ?he 
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Finally, the quasi-steady and radiation-linearization 
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2-d-i Model Equatinnc} 


at^ither o^tttMdH l” *7^ ^ ^ ^ 
in these annularM^d ^*^i»" ~ ' 
temperature field, but without feedback to th. ,h ? ^ 

For small Rayleigh numbers and an infinitely^ Lng a^ulu^^Uh 


ORIGWOt 

ps: p00* ? . QbW-lTt 


2 


Thui 1 f:r lta “° nal accele "-ation [Arpaci and Larsan 1984] 
.the conduction heat transfer from surface i to im 
Vi * 1 or 2 ) is given by: 


Q ki - 


T,- - T ; , 


where resistances are 




> R[ 




T 1m - T- 


, R| 




(3.2,3) 


and the effective liquid conductivity is Similarly 


°h 2 ~ h^A 2 ( T 2 ~ Tf ) 

Energy balances on surfaces 1 and 2 yield, respectively 


Q c = Qm + Q. 


1 ' Q rl “ °k2 + Q h2 


(3.4,5) 


(3.6) 


(3.7,8) 


where Q rl is given by (2.15), and (2 16) or 171 «.u . 

bound"- 0 " ?h e h8at 3dded reSUUS in the * the 6 

boundaries through H s& , the latent heat: 


Q kl ” Q ks “ ~ 

at 


(3.9) 


°k2 + Q ks “ + 2 «^ h s £ r 2rn~ 2m 

dt 

— - - 2 Adiaba ti c ...Inner Wall ( Case a ) 


(3.10) 


For the adiabatic inner wall, with R, = r = M fk. u 
balances and rate equations combine as follows: ’ ° V6 ener9y 


f A U C 1 

( 1 - A5 1 )j- A - = i - n( e. - ) 


(3.11 ) 
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(3.12) 


( r + a s 2 )— 2- = o + nc e - e ) 
dr i 2 7 




N( $2 “ = 0 


where 5 ± and S 2 are the scaled liquid-gap thickn 


esses 
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1 ” r. /r . 


r l ' r 2 


*2 = 


r 2m r 2 t 0m /r -i — r 


r l ' r 2 


the scaled excess temperatures (i = 1,2) are, 

6 i - < T X - T m > /T r 
the logarithmic terms are, 

v = &a{l/( )}/A, u » &i{1 + A5 2 /D/A 

the radiation number is. 


N '' TT m KrlCri r 2 )/k e&" h rl^ r l ~ r 2 )/k e£ (3 

where r - t/t r , T = r 2/ r jL , A = 1 - r, and where the reference 
temperature and time magnitudes, T r and t r , are determined as 


T r = 


PrCr, - r 0 ) 


V = 


( r 1 ~ r ■> )P/)H. 


J 1 , . q r “ q s for case a ■ Typical parameter values are shown in 
3 6 3 1 ' Xt is noted that the time constant represents the 

ratio of the mass of PCM times the latent heat, divided by the 
source flux applied at the arithmetic average radius. The 
initial conditions for case a are that all ^(0) = 5 .( 0 ) = 0 . 


4 


ORIGINAL PAGE IS 
OF POCtf QUALITY 


3j3 - Xaayea ^g-Xnnar ui aii m... tg 


For the convective boundary <-or„~r*- 

R k 2 ■ o CT 2m . T ) and T ar r « *" •* the inner W a U . ulth 
equations combine as follo us: m f ’ the ener9 >' fiances and rate 


(l - A 5. ) 


2!] 
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s s - m i + $ 2 ) 


(3.23) 


(3.24 ) 


(3.25) 


V " 9 r 

N< ° 2 ' ) * * ■ - #rc i + # 2 ) 

«here the additionai losarithmic term ls 
w = ^{( 1-A5 1 )/r>/A 
the effective conductivity ratio is y u , 

15 fi = h f (r l ~ r 2> /k e4> and where q is oh>? ^ th * Bi0t number 

ei. Here q r ls obtained from (3.21). 
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4. RESULTS 


l^l-^^_batic_JjTn er Mall ( Case «) 


From (3.11) and (3.12) it is seen t har - 

of the radiation number: H heir SUm 1S inde Pendent 


(1 - AS* )2i 

1 dr 


+ ( r + A5 0 ) 


d5- 

— < 

dr 


= l 


(4.1 ) 


which has the exact integral 


*1 + 


rS 2 - J [»1 2 - «**) - 7 


(4.2) 


In particular for 5 + ^ — i i_ 

melted, (42) yields the 2 m ~ WhSn ^ Solid is completely 
• ; yields the maximum r, or scaled melting time 


T m = 


i-LE = Ll. + r 2 

2 2r, 


(4.3) 


charging or ng't'he PCM "in o^tai^d a T 
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y o vo.14 j can be combined as 


d ( ) N5. 


d( i+nt ) ' 1 + Nr 1 

which. With (4.2), has the solutions 


(4.5) 
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N5 1 = 


2 + NT 

NT 

NT 


, Ntf 0 = — 


1 + NT 

2 2 

1 + NT 


(4.6,7) 


N<?^ = NT 


2 + NT 1 + ( 1+NT) 2 
4 ( 1 + NT') 3 


, N $2 = 


2 + NT 
4 


NT 

1 + NT 


(4.8,9) 


These solutions are shown in figures 4-1 and 4-2, where the 
double surface phase-change phenomenon is clearly evident. It is 
seen in (4.6,7), and in figure 4-1, that S grows linearly with T 
for small NT(i.e., initially), whereas $ 2 grows quadratically ; 
moreover, the growth of 5 2 depends on the radiation number, N, as 
seen by the ratio (4.10) of (4.7) to (4.6) 


SZ - 


NT 

2 + NT 


2a. 


( NT) 2 

1 + ( 1+N-T) 2 


(4.10,11 ) 


where S 2 >0 for N >0, and 5 2 > 5 ^ for N > co . That 

is, for small radiation number the phase change will occur 
Primarily from side 1, whereas for a large radiation number it 
will occur equally from both sides; for N = 1 and the end of the 
charging process where in the present case T = 1, the ratio of 
Phase growths in (4.10) is 1:3. However, depending on the 
surface emissivities , N and the radiation effect can be quite 
small, as indicated in Table 3-1. 

In the present formulation, wall temperatures are determined 
rather than specified, as in figure 4-2. The gradual increase in 
the wall temperatures (absence of thermal boundary layers) 
validates the quasi-steady approximation, and with T = 26 6 K 
the small changes of Tj> and T 2 from T m validates linear ization ’of 
the radiation term. In (4.8,9,11), and in figure 4-2, it is seen 
that initially increases linearly with T, but that 
increases cubically with T. 2 

In the case of an annulus with finite radii ( r /r = r y* l) 
equations (3.19) to (3.22) were solved by numerical integration, 
with ^results as shown in figure 4-3, to the same scale as 
previously, but for T = 0.5. Here the phase front movements were 
calculated for the series of radiation numbers shown, up to the 
maximum time +, = (l + D/2 = 3/4. There appears to be no 
dependence of NS 2 on N, and only slight dependence for N^ ; 
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however, the curves are shifted somewhat, relative to the "planar 
annulus" result in figure 4-1. Without the N-scaling, the phase 
fronts appear as in figure 4-4 for the radius ratio, 0.5. 

4.2 Convective Inner Wall (Case b) 


With the convection heat sink at the inner wall, the possibility 
exists that q g may not be sufficient to cause any melting at all. 
This is seen in (3.23) where the initial slope must be positive 
for to grow; thus, g s must be greater than the minimum 


N + KyWfv 

9 so = *rxi + <? 20 ) = fir 0- 


n + k:/w n + nr 


(4.12) 


where (? 2 o obtained from (3.25) with 0 
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= 0 : 


^20 “ ^2^ 0 ^ ~ 


- fiT 


N + K/w n + fir 


(4.13) 


where w Q - lim(w) as 3^ — > 0. In dimensional variables with 

^rl = crT m ^rl ’ limiting minimum heat flux for phase change to 

occur is 



+ II 

h f r 2 h rl r 1 + k s /^.{r 1 /r 2 ) 


which is the initial overall temperature difference over the 
total thermal resistance. 

Similarly, there is a critical heat flux for which the PCM just 
becomes 100% melted and at steady state with zero slope at 
<5^ - 1. This is given by (3.23) as: 

9 sl “ ' er (4.16) 

which, in dimensional terms, simply is 

q sl = h f (T m " T f )r 2 /r l (4.17) 

Finally, it may be asked if there is a q g above which melting 
also occurs at the inner wall. This limit is given by (3.24) and 
(3.25) with $2 = 0 : 
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= fiT 
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However, for small N as in Table 3-1, this flux limit is quite 
high so that (unlike case a) two-surface melting does not 
usually occur with the convective boundary condition and 
practical flux levels. 


To summarize these conditions we have case bO ( 9 S <_ 9 s0 > which 
results in no phase change, bl ( 9 s0 < 9 S < 9 sl ) which results in 
a steady state with only partial melting, b2 (9 S ^ 9 S !-9 S 2) 
which results in 100% melting in finite time from the outer 
surface only, and b3 ( g s > 9 s2 ) which results in melting from 
both surfaces. 


The numerical integration of (3.23) to (3.25) for case bl is 
shown in figure 4-5, where the scaled heat flux was taken as 


g s = *g sl + (i ~ *)9 s0 . 0 <-1 

The partially melted, steady state condition is evident, as well 
as the effect of the surface emissivity. For comparison, 
computations were also made with the full nonlinear radiation 
term (2.16), but the results differed less than 3% from those in 
figure 4-5. Partially melted, and unmelted, conditions have also 
been found experimentally [Strumpf and Coombs 1990] . 


This steady state condition may also be seen analytically for the 
parallel plate limit, A — > 0, where the system equations combine 

as follows 


dx £ ( N-x )( g^x+N ) + KN(l + x) g (4.20) 

Nd r + N ( N-x ) [1 + £(l + x)/N] + K(l + x) S 

where x = NS,; this is a quadratic form which, with (4.19), has 
the steady state solution 


<K N + K + ft) 

1 ,ss + N + K + ( 1 — 0)( K-l )N 


(4.21 ) 


For case b2 the phase-change front grows rapidly to complete 
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melting, as shown figure 4-6, where the applied flux was taken as 
9 S = *8si. Also shown is the effect of radiation, which may be 

M 0t " e : m ° delS “ here the effect of radiation was not 
included [Kerslake and Ibrahim 1990] 

The discharge mode, in the present context, utilizes the same 

r :: r u t- (3 -«> 0.2s), bU t w ith 9s = 0 . Th e r .s U i t 

s the decay of ^ to 0, from the melt condition when the applied 
flux is turned off, as shown in figure 4-7 for several case b2 
conditions. 


CONCLUSION 

The Phase-change behavior of a high-temperature salt enclosed in 
an annulus with specified outer-wall heat flux has been analyzed 
including the effects of radiation within the enclosure and two 
inner-wall boundary conditions. 

For the adiabatic inner wall condition, radiation resulted in the 

melting of solid PCM from two surfaces, and the time to complete 
melting is a fixed quantity. 

For the convective inner wall condition, melting was found to 

occur from the outer wall, only (or not at all), for practical 

heat fluxes; partially melted or fully melted conditions were 

found depending on the flux level; the effect of radiation is to 

by pass the melting process and, thus, increase the time for 
complete melting. 

For both boundary conditions, the imposed flux condition resulted 
gra ual changes in the wall temperatures, thus validating the 
quasi steady models. This contrasts with the traditional Stefan 
Problem with imposed rapid changes of the wall temperatures. 

Radiative exchange between the two walls was found to be about 
twice as high for the LiF medium as for vacuum; however, wall 
temperatures differed only slightly from the melting temperature 
during the phase-change processes, and the overall effect of 
radiation was about ten percent of the total heat transfer rate 
for practical heat fluxes and annuli dimensions. This results 
from the assumption of continuous contact of the liquid PCM with 
t e canister walls; further analysis is required to evaluate the 

effect of voids within the canister, where radiation would be 
more important. 
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Table 3-1. Typical Parameter Magnitudes 
Parameter Flight Exp.* Space Station* * 


Outer Radius, r^ [cm] 

3.51 

2.11 


Inner Radius, r 2 [cm] 

1 .95 

1 .19 


PCM, Melting Temp., T m [K] 

LiF, 1120 

LiF-CaF 

2 , 1040 

Latent Heat, H s ^ [J/g] 

1037 

816 


Liquid Density, pg [g/cm ] 

1 .79 

2.19 


Solid Density, P s [g/cm 3 ] 

2.33 

2.59 


Liquid Conductivity, k^ [W/cmK] 

0.037 

0.017 


Solid Conductivity, k s [W/cmK] 

0.060 

0.038 


Radiation Cond., k r [W/cmK] 

0.020 

0.018 


Radiation Nr., N ( £ . =0.1/0. 3) 

Ml j V 

0.102/0.365 0.069/0 

.242 

Heat Source, q s [W/cm ] (avg./peak) 

0.921/ 

0.650/1 

.05 

Fluid Film Coef f . , h 2 [W/cm 2 K] 

0 

0 .028 


Biot Number , & 

0 

0.81 


Fluid Temp., T^ [K] (min max) 

— 

800 

1000 

Temperature Scale, T r [K] 

26.6 

240 

40 

Time Scale, t r [minutes] 

52.4 

3.0 

18 

Heat Flux Limits [W/cm 2 ] 
q s0 ^ t m,v =0 ’ 1/0 ' 3 ^ 

q s2 < c -m,v=°- 1/0 - 3 > 


1 .75/2.95 
3.79 

45.5/15.7 : 

.29/ .49 
.63 

7 . 6/2.6 


* Namkoong 1989; Williams 1988. 

# Strumpf and Coombs 1990; Kerslake and Ibrahim 1990. 
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1* Spectral Refractive Index for Lithium Flouride, 
and Blackbody Emissive Power at 1120 K, 
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Fig. 4-4. Growth of Phase-Change Fronts for Annulus, 
Adiabatic Inner Wall. 
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Fig. 4-5. 


°f Phase-Change Front for Partially Melted 
Condition, Convective Inner Wall. (Annulus 
parameters per Space Station, Table 3-1.) 
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e ' 4 6. Growth of Phase-Change Front for 100* Melted 
Condition, Convective Inner Wall. (Annulus 
parameters per Space Station, Table 3-1.) 
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1 - INTRODUCTION 


to 


understanding^!? ^hase-chang^ther *°1 * S Central to the 

development of Solar Dynami? f 1 ene ^ 9 / storage for 

Shuttle flight expe^imenfo k Pa ! e Power • Therefore, Space 
1989 a , 1989 b], and numer ica f*** 2 ®®" planned fNamkoong 
determine two and three dimeriitZTtf b ®! n perfo ™*d 
transition process [Kerslake Ulo J f !® CtS of the phase 
analyses and experiments utilirl ’ Wlchn f r 1988 ]. These 
containing the phase chanoo la f f n annular canister 
heat flux is Impressed on t h r te : ial (PCM) * where ■ solar 
( here radius r , Ttomel t t PcTld "^ 1 ° f th * «nister 
inner wall (here radius r-O tofr-II* Conve f tlon at the 
(in the experiments the solar fil !” 6 °[ solldlf y the PCM 
simulated, and the comLI?« T lux . ls electrically 

radiation). This annular c^niste^ae ^?' * ^ conduct ion and 
current design for the Solar nun geometry is also the 
Proposed Space Station ° f th ® 

temperature silts I^ch^^lithii^f ? th ? PCM ’ S ar e high- 
melting temperature of nil i r f ^ oride ( L i p ) with a 
temperatures radiative tranJ! \ 1588 F) ’. At these 
the overall heat transfer Drill 1 1S a . S1 ? ni f icant part of 
In the previous two-dimensiona )' SS ®2 ^ lbh i n tbe canisters. 
1990] radiation was - [K ®:? lak ? and Ibrahim 

dimensional model [Wichnar i ,2 \ n the Three- 

radiation is obscured by t ha I a1 ,* 1 ? 88 ^ the effect of 
computations; in the fliahi- ^ omp l? xl ty of the numerical 
radiation will occur natula l l expar ifner »ts [Namkoong 1989] 
cannot be deterged dxnectl^as 0 ^^* r ? distlve 
temperatures will be measured- !i ^ canister surface 
radiation in vacuum ua C * > although convection and 

’*>9) the combined ^rad!ai?on Stl 3 at ! d Gently [Yucel.t el., 
generally treated in the technical^?? Ch f n9e effec t is not 
, a nee d for basic models Ind III ^ terat “ re - Thus, there 
this effect. To this eldVll 2- " S whlch illustrate 
Presented. d 3 on e-dimensional analysis is here 


actual U heat H f luxes S may S be nv f r are finite in length 

ccrform, there are a number 8nd anally 


why 


a one- 


and 

ci menIioIal ’aDD ere - 3re - 3 number of reasons 
^ PP L°i lmation is useful: 

which i 1 1 ust rate 1 the s °l utions ar e possible 
interactions ; 3te 3S1C pba nomenological 

numerical moSels'ba y^be a9alnst >- h ich complex 
boundary condtuons- COmpa ’' ed identical 

° estimation of experimental 

variations; surface temperature-time 


b) 
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dimensional ^in^th^radial - pre ^? m * nant ly one- 

«th larger 

convent ional °ex change betwee^diff enclosure Allows 

, . ? _ i n ^ erver| i ng medium, LiF hse USe Sur ^ aces » except that 
d^fer from those of v^cuim’ ?k SPactral Properties which 
application requires spectl'i T h f refoir e, accurate 
Williams [1988] a two-band nt ? 9 ra ti°n . Following 
recent measurements of LiF optici^ mati ° n is made uti ^i?ing 
Hunter 1985]. it is sho^n ?n c^- Pr ° Perties t Pa Hk and 
temperature, T m = 1121 K th " ^ c . tl ?! 1 2 tha t at the melting 

For e ^h° CCUr f, in the tra nsparent pa^of^ radiation emissive 
or the smaller part at lnnoo Part of spectrum for LiF 

approximation is used which !dds a »Q en8ths- a * th ick-gas ; 

thermai conductivity [Williams^ MeJ ° r l ° the appare "t 

is assumed to follow the hLI^d metal cani ster surfaces 
Howell 1981, P . 112] . with t^ en ? r f latio " Siegel and 
this leads to a spectral ] v inf s Peotral properties of LiF 
function which is the riiio of^hl^ emissivit >' fractional 
surface contacting LiF to the\Jw ? v ?rage emissivity of the 
vacuum at the same tempe?a[ure °t the »urf« e in 

energy exchange model is then deve l two ~ su r face radiant 
linearized, results in a fart J ? ° Ped ? hich ’ when ■ 
radiation heat transfer coefficient 5 ' 47 10 the ec ^ uiva lent 
usual factor of 4 , a 37 ^ i lent ’ as compared to the 

increase over radiation exchange in 

3, wherl the PCM^S i ij* ^ fi y M °t e 1 ® consldared In Section 
flax is suddenly applied at !hIV " Peratura - T "" and a heat 
“file the inner surface r.l?J canister outer surface! 

Problem is one in which the f S abla batic. The usual Stefan 
suddenly changed [e g Jao ^ P * rature at a boundary is 
° T . Solomon 1981] , relultino Pyusa 1989 * Burmeister f 9 83 
as the error functo r ? b f^dary layer growing in ’ 

?Sx L £S3.ir 43 * 

'emaiI 1Cantly and suc h that the pha^ layers do not develop 
mams axially uniform T h^ • h f' chan9e boundary 

* 985 ^ 1700, for rectangular enclo< 5 valld for Rayleigh numbers 
^! 5 > P- 401 ], but has also Z21 loaures Tlncropera and DeWit? 

and an° h lar9er Ra/leigh numbers with nSt h ated fr eXPerimentaI ^ 
tn e f?- a f Pect ratio of 4.5 fz L ^ t he3t flux boundary 

1 9 t experiments the Rayleigh n numh Jan 198 ?^ ’ since in 

ayieign number remains small 
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of C natnr = t° n ' Un ?® r n, f cr °9 r a^ity conditions the'mt 
convection tffeltl " A f suppressed - although some ~ 

Arnold er a !! !?90] may PTeSent ^ a ^ 9 lu and Lam 1987 

n sumptions regarding the second boundary condition Cat 
r^) are tied in with the radiation model (i.e., withou? 

uinul^H aner9y e *change between the two walls, the PCM 

would Simply melt from the outer surface and the effect 

front ■a°rr?v~- r tK t ^ \"" er SLjrfa ~ urt11 theph^-chanse 
l°2l n ! S there). Here, an adiabatic condition is 
chosen as representing the best that can be done C as far a* 

results ‘ '"other > c yieldi "9 the analytical 

esults, other conditions, such as convective or thermal- 

to a large W exJent han9e ^ PreSent results somewhat, but not 

inner T w!l [ 6S ti tS °f th f analysis are that, with an adiabatic 
“ !1, *; he rate of melting of the PCM is independent of 
adiative and conductive effects within the annulus but 
depends only on the heat flux at the outer wall and the 
olume and heat of fusion of the PCM. The effect of 
radiation and liquid conduction is to redistribute energy 
and cause melting of the solid PCM at both its outer and 

Process it^as id . con ^ uc tivity did not enter into this 
f ocess It was also found that the wall temperatures 

"'idaUno Sm °°) hl C and Slowl> - durins the Pel tin^r thus 
^a.idating initial assumptions. 
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2. TWO-BAND RADIATION MODEL 


To determine the effect of radiation within the LiF 
enclosure, it is necessary to know the optical properties of LiF, 
which differ considerably from vacuum. In particular, if the 
specular optical constants are known ( refractive index, n, and 
absorptive index, or extinction coefficient, k ), other radiative 
properties can be obtained. A summary of recent measurements 
[Palik and Hunter 1985] is shown in Fig. 2-1. It is seen that 
there is a transparent region of the spectrum near 1 /an (where 
k — > 0), which is the region of interest for emission near the 
melting temperature of 1121 K. The refraction index appears 
relatively constant in this region, but decreases somewhat with 
increasing wave length. This is shown in Fig. 2-2 where the data 
of Palik and Hunter are replotted on a linear scale. Also 
plotted is the Planck spectral emissive power for 1121 K, which 
shows that emissions whithin the canister lie primarily between 1 
/.on and 6 /.mi . The corresponding extinction coefficient data are 
shown in Fig. 2-3. Thus for most of the relevant spectrum, the 
single c .rystal LiF is transparent to radiation. This is also 
shown by previous data [Amr . Inst. Physics Handbook 1957] , as 
reported by Williams [1988] . 

2 . 1 Blackbody Emission 

For the transparent part of the spectrum, "K <_ , the 

blackbody emissive power is given by 



where c^ and c 2 are the radiation constants (based on c 0 , the 
speed of light in vacuum): 

c 1 = 2-7thc 0 2 = 3.7415 x 10 -16 Jm 2 /s 
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c 2 = hcQ/k = 14,388 iflnK 


The index of refraction for LiF is the spectral quantity shown in 
Fig. 2-2, which is here approximated as 


n = 1.39 - 0.00271X 2 * a - b\ 2 
To evaluate (2.1) the change of variable is made. 


£ = c 2 /n>,T 


so that 

d£ = -£ 


n + Xn ' 


rVX 


dX 


( 2 . 2 ) 


(2.3) 


(2.4) 


where n' = dn/dT = -2bX. Substitution of (2.3) and (2.4) into 
(2.1) then yields the expression 


' b i s 


t 4 15 

•7T -Jf 

■n 


.CD 


£-d£ 


(2.5) 


+ >,n' e' - 1 


where £ a = = C 3 /T = 1995/T (K) . 


An example of the numerical evaluation (2.5) is shown in 
Fig. 2-4 for > >a = 5.5 £tm , where the upper limit of £ = 16 is well 
whithin the validity, of (2.2). This may also be written in terms 
of an effective index 


2 4 

e = n cr T 
b ,a a 


( 2 . 6 ) 


The effective index, n a 2 , obtained from the integral in (2.5) 
with = 5.5 A^m, has a temperature dependence because of the 
spectral dependence of n( T ) , as shown in Fig. 2—5. This result 
is similar to that previously obtained [Williams 1988] , except 
that previously the effect of n’ (a 10% effect) was not accounted 
for and the blackbody fractional function was included in the 
effective index. 

2.2 Emissivitv and Absorptivity 

The assumption is made in this analysis that the metal 
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surfaces bounding the LiF are diffuse. For spectral emissivity, 
the emissive power of the surface is given by 


n. 


e a = 




so that with the above change of variables 


e. 


a T 


4 15 


s ft 




n + Xn * 


£ 3 d£ 

^7 


(2.7) 


For metals, a recommended formula for the normal spectral 
emissivity is [Sparrow and Cess 1978] 


'Xn 


= 36 . 5\|r/X ’ - 464 r/X 


( 2 . 8 ) 


where r is the electrical resistivity in Ohm-cm , and X is the 
wavelength in /.on. At "X * 2.5 (.m (maximum radiation at 1121K) and 
the rather large value r = 10 3 cm, the second term is only 
2.5% of the first term and therefore neglected. The result in 
vacuum of spectrally integrating (2.8) is total normal emissivity 


£ n,v = 0 - 576 J rT s ' (2.9) 

where T s is the surface temperature in degrees Kelvin. Let r be 
proportional to the absolute temperature, then r = ’TqT s /T 0 , where 
r Q is the resistivity at temperature T 0 , and where from (2.9) the 
total normal vacuum emissivity at this temperature is 
Substitution of these variables into (2.8) results in the normal 
(or hemispherical) emissivity ratio 



6 0 , v 


36.5 

0.576T 0 


•\T 


v 1 



where c 4 is the nondimensional number 


( 2 . 10 ) 


c 4 = 36 .5/0. 576 \| c 2 ' = 0.529 
Now, substitution of (2.10) into (2.7) results in 


of po or quality 


1, f T I " 


( 2.11 ) 


3 a T Jvvl £ -o,v wv 


where the emissivity fractional function is defined by 


f t.( T ) 



,nOD 


n 


3 


n + \n* 

J ^a 



£ 3 d£ 

^7 


( 2 . 12 ) 


This function is shown in Fig. 2-6 for = 5.5 A<m , where it is 
seen that the effective total emissivity is roughly twice that of 
the value in vacuum: 


- Vv < T S /T 0> f £a<V 


(2.13) 


A completely similar procedure yields the effective absorptivity 


“a< W * 

where T g is the 
2-3 Thi ck Gas 


^0 i v 


T_T 
s e 


environment 

Approximation 


‘ f Ca<V 

(other surface) 


temperature . 


(2.14) 


It would appear from Figs. 2-1 or 2-3 that k is sufficiently 
small for the wavelengths of interest that the medium can be 
considered perfectly transparent. there is, however, a region 
at the longer wavelengths where absorption occurs in the medium, 
as can be shown by evaluating the absorption coefficient [Siege-1 
and Howell 1981, p 427] from the data of Palik and Hunter [1985]: 


a. 


A 


4nk 

X 


as shown in Fig. 2-7. These values may also be read as the 
optical thickness for a 1 cm layer of LiF (or as half the optical 
thickness for a 2 cm layer). The results are in general 
agreement with previous measurements [Amr . Inst. Phys. Handbook 
1^57], thus, the assumption of an optically thick gas above >, = 6 
0r 7 /.mi is supported by the more recent measurements. Therefore, 
following Uilliams [1988] , the two-band approximation can be made 
that the LiF is transparent to radiation below wavelengths of 


GRiGSftAL PAGE IS 

OF P00K QUALITY 


approximately S.S m. and acts as a mtl. 

^ :sr - i t h - surfaee prop 

2 4 R^d i&rjt Ener^ xJ^bapgg- 

For each surface, i. it an enclosure we have the net 
radiative flux. 


(2.15) 


q = c l. " q : 
i i 1 


- J if f lise opaque surfaces, 
and th . radiosity definition for dtffus . 


(2.16) 


q = C B + P.<*. 
q i i i 1 1 


4 is th e blackbody radiosity (2.16) we 

where B- = * T i ^ - between (2.15) and 

Eliminating the irradiance, Q, . ° 

have 

a - u e. - « < *U /P , 


or 


+ 

c - q. 
i i 


(2.17) 


0 = 


l 

• q - q-A- , the modified 

where the total heat transfer rate is^j sur face 

kiarkbody radiosity is c i 11 1 


blackbody 
resistance is 

1 " o: . 


( 2 . 18 ) 


R = a 

i O'. A 

i i 


Additionally, the follow 

medium , 


coolies in the transparent 
ing relation applies 
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N 


= > 


J=1 


a f .(q. 

i i J 1 


N 

+ vn 

q ) = ,> 

V ^ 

j=i 


+ + 

q - q 
1 J 


(2.19) 


1J 


so that conventional radiation ^“®^ 0 ®® l ^d°with the surface 

the modified blackbody surface orm ana lytical 

resistance given by (2-«>. * ^ daper , dE on the -other-surface 

is. Of course. not po!! J (2 .14)1 however . an adequate iterative 

temperatures, as s ow weighted average for T e - 

solution can be obtained wxth a weig 

^ * i 


FO r a two-surface enciosure, (2.17) to (2.18, Xi.U 


C l' C 2 


Q = 
1 


- Q 


1-0 


A 0. 

1 1 


A F 
1 12 


1-0 


A O 
2 2 


( 2 . 20 ) 


where , with (2.13) and (2.14). 
g(T ) 


l 4 

— a T 

i g( T ) 1 


c = 


g(T ) a 

2 T 

r rr >7 T 

2 g ( t ) 2 


g(T) - f. (D'-f^ 
ba 


( 2.21 , 22 ) 


(2.23) 


a linearization of (2.20), C x and C 2 may be expande 
JoTii- Lout the melting temperature. V 


C * C m + (T 1 T m )C i,l"> + (T 2 T,1,) i,2m 


dith the result that 


r _ c j- (4 + 2T g V 9 ) OT (T T > 
C l C 2 m m m "i 1 

or. with (Tg */g) m - 0.237 from Fig. 


( 2 .24 ) 
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(2.25) 


c - c - 5.47 OT (T “ T ) 

1 2 ml ^ 

Thus, the radiant heat transfer rate between two surfaces may be 
approximated as 

. _ x (2.26) 

0 = - 0 = h A( T " T ) 

l 2 rl 1 1 2 

where the effective radiation heat transfer coefficient is 


5.47 cr T 


HL 


(2.27) 


rl 


1 - 0 ! 1 A 
-i- + + 


1 - 


Hz 


Cta 


F 12 A 2 


ck 


• a. x 07 ^ Hoes not involve differences in 
?ince the denominator of (2.2 ) evaluate a, and ou at 

Konrot i vities it is considered adequate to evaluate a* ^ 2 
absorpti vi utilized subsequently m 

the two-surface Phase chanse 
problem. That is, for Oj = «2 * , and T m 


rl 


5 .47 a T 


m 


43.7 tJ/nv K 


■'m 



1 r 


[l + r 1 /r 2 J - 1 


- + r l /r 2.l 


(2.28) 
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3. BASIC PCM CHARGING MODE 


The LiF PCM contained by the annular canister in the flight 
experiment, as well as in applications with annular geometry, is 
to cycle between the charging (heat addition) and discharging 
(heat removal) modes. In the first mode heat is added to 
(impressed upon) the large diameter, external wall to melt the 
PCM while some heat may or may not be removed from the 
small-diameter, internal wall; in the second mode, heat is 
removed from the interior wall to solidify the PCM while the 
outside wall is, more or less, insulated. 


Here is considered the basic charging mode where at time zero the 
external heat source is Q.(W/m) and the PCM is 100% solid and at 
uniform temperature, T m , and where the inside wall remains 
adiabatic throughout the charging process. Furthermore, to 
obtain a basic (or fundamental) solution to the charging problem, 
an infinitely long cylinder is assumed together with axissymetric 

heat input. 

The above assumptions result in a radially one dimensional 
time-dependent system, as illustrated in Fig. 3-1. because of the 
transparent part of the LiF spectrum there is radiation from the 
outer wall (surface 1) to the inner wall (surface 2), where the 
latter acts as a refractory or reradiation surface, such that t 
energy source, Q_ , causes phase change at both melt-surfaces, lm 
and 2m, as shown. Additionally, axial gravitational acceleration 
is assumed such that any void is removed to the end region, an 
such that liquid remains continuously in contact with the two 
walls and the phase-change boundaries. 

Although the above stated problem is quite restrictive from a 
practical view, it nevertheless results in closed-form solutions 
to which more general numerical solutions maybe c ° m P are f|; 
Additionally it enhances our basic understanding of the PCM 
charging process when radiation is included. 


3 . 1 Model Equations 


The model considered is for a solid region surrounded by liquid 

as indicated in Fig. 3-1. In these 


at both radii r and 


2m’ 


annular liquid gaps convection is 
transport equation 


described by the general energy 


1 
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dl 

St 


(3.1) 


ST ST o , - , n 

■ + u — + v — = aV^T (.3.1; 

Sx Sr 

where u,x are in the axial direction and v,r are in the radial 
direction. However, under the parallel flow assumption for long 
cylinders, u = u(t,r), v = 0, and T = T(t,r). This means that 
the fluid transport terms on the left hand side of (3.1) are zero 
and the temperature field is determined entirely by conduction 
(including the conductive effect of the thick-gas radiation 
absorption), even though there will be convective flow driven by 
the Boussinesque temperature term [Arpaci and Larsen 1984] : 


v S(ru) Su 

g£( T - T 0 ) 

r Sr St 


(3.2) 


Thus, the long-cylinder, parallel flow assumption essentially 
uncouples the momentum and energy equations such that the 
temperature field is determined by conduction, independent of the 
velocity field, and then can be used to determine the velocity 
field as in (3.2). 

Because the movement of the phase fronts is much slower than the 
liquid response the problem may be regarded as quasi-steady with 
known solutions utilized for the temperature distributions. 

Thus, the conduction heat transfer from surface i to im for (i = 

1 or 2) is given by: 


T: “ T, 


(3.3) 


where resistances are 


&nt r i /r « m ) 


. R, 


£n(r 0 /r^>) 


(3.4,5) 


and the effective liquid conductivity is k e £. An energy balance 
on surface 1 yields 


°s = °kl + °r: 


(3.6) 


where Q rl is given by (2.26), and for the adiabatic surface 2, 
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0.7) 


CL 1 =0, 
rl u k 2 


0n a flux basis (2 26) f-; ^ 

’’ C - 3> - Snd (3 ' 4) *" d <3.6, cosine as 


h rl (T l ~ T ) + ^- (T 1 ~ T^) 

» r. 


(3.8 


m J 

r !&Xri/r lm ) " q s 

Snd ^ 2 - 2b '>- (3.3). (3.5) and (3 7 ) r „ . . 
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It is noted that the time constant represents the ratio of the 
mass of PCM times the latent heat, divided by the input flux ^ 
applied at the arithmetic average radius. The initial conditions 

are that all 0^0) = 5^0) = 0. 


3.3 Solutions 

The system, (3.19) to (3.22) does not possess a closed-form 
general analytical solution (primarily, because of the curvature 
effect represented by the logarithmic terms). However, by _ 
combining (3.19) with (3.21), and (3.20) with (3.22) we obtain 
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Hence, with the definition of the time constant, the time for 
charging or melting the PCM in the annulus is obtained as 


* 


” ( ri 2 - r 2 2 >^. H sft 
2wr l q s 


(3.28) 


that is, the melting time is precisely the phase-change heat 
content divided by the total rate of heat input, independent of 
radiative and conductive properties! This is perhaps a 
surprising result, but quite reasonable: since in the basic 

charging model the phase-change boundaries are the only sinks for 
the energy source, 2rrr 1 q s , the net result must be the melting of 
the PCM, with conductive and radiative effects redistributing the 
energy input and determining the relative locations of phase 
boundaries . 


The above behavior can be seen analytically for the limiting 

condition € > 0 (F > 1) which is the parallel plate limit. 

In this limit, the system (3.19) to (3.22) becomes 
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such that substitution of (3.33) and (3.34) into (3.29) and 
(3.30) yields equations for the phase fronts, only: 
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One solution is obtained from the sum of (3.35) and (3.36), or 
from (3.26) with 6 = 0: 

Si + S 2 = T ‘ (3.37) 

Eliminating S 2 between (3.37) and (3.35) yields the Euler 
equation in 1 + NT, 
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Thus, the solution to (3.37) and (3.38) is 
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It is seen in these equations, and in Fig. 3-2 that 5^ grows 
linearly with T for small NT (i.e., initially), whereas S 2 grows 
quadratical ly ; moreover, the growth of 5 2 depends on the 
radiation Biot number, N, as seen by the ratio of (3.40) to 
( 3.39 ) 


£2 NT 

5^ 2 + NT 


(3.41 ) 


where 5 2 > 0 for N > 0, and. S 2 > 5 ^ for N > to . That 

is, for small radiation Biot number the phase change will occur 
primarily from side 1, whereas for a large radiation number it 
will occur equally from both sides; for N = 1 and the end of the 
charging process where in the present case T = 1, the ratio of 
phase growths in (3.41) is 1:3, or, 75% melted from the outside 
surface and 25% from the inside surface. This result contrasts 
with other recent models where the effect of radiation was not 
included [Kerslake and Ibrahim 1990]. However, regardless of the 
value of N, the PCM is 100% melted when T = 1 (for this parallel 
plate limit of the annulus); this is a consequence of the imposed 


heat flux boundary condition arising from actual or electrically 
simulated solar irradiation, which differs from the Stephan-type 
problem where the bounding temperatures are specified, and from 
the assumed adiabatic inner radius. 

In the present formulation, wall temperatures are determined 
rather than specified. Thus, substitution of (3.39) and (3.40) 
into (3.33) and (3.34) yields the temperature functions 
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and the ratio 
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Here, and in Fig. 3-3, it is seen that initially 0^ increases 
linearly with T, but that <? 2 increases cubically with -T'. 

In the case of an annulus with finite radius ratio, r 2 /r^ = ^ » 
equations (3.19) to (3.22) were solved by numerical integration, 
with results as shown in Fig. 3-4 with the same scale as 
previously, but for T = 0.5. Here the phase front movements were 
calculated for the series of radiation numbers shown, up to the 
maximum time T' m = (l + D/2 = 3/4. There appears to be no 
dependence of N5 2 on N, and only slight dependence for NS^; 
however, the curves are shifted somewhat, relative to the "planar 
annulus" result in Fig. 3-2. Without the N-scaling, the phase 
fronts appear as in Fig. 3-5 for the radius ratio, 0.5. 


The results of the present analytical one-dimensional model may 
be compared with those of the numerical three-dimensional model , 
NORVEX . Such a comparison is shown in Fig. 3-6 using the 
following data: 


3.51 cm, r 2 = 1.95 cm, T m 


1121 K, H = 1037 J/g, 

Sjf 
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= 1.79 g/cm 3 , j'-' s = 2.33 g/cnw, 
k^ = 0.037 W/cmK, k c = 0.06 W/cmK, <k e =0. 037+0. 017=0. 054) 

= 2.45 J/gK, c 5 = 2.45 J/gK 

fl = 0.3, N = 0.15, h rl = 0.00524 W/cm^K, q £ = 0.921 W/cm^ 

where the canister length was 6.6 cm, and where the solid 
unmarked curves are from the analytical model and the marked 
curves are from the NORVEX calculations. It is seen that 
analytically a smooth increase in outer wall temperature is 
predicted, whereas the numerical results show no increase above 
the initial melting temperature in the first 40 minutes of 
heating. Also, the analytical melting rate- is somewhat higher 
than the numerically predicted rate. 

These differences are not easily ascribed to three-dimensional 
versus one-dimensional effects as only weak axial gradients 
occured in the numerical calculations, where the canister aspect 
ratio was 4.2. Nor do the differences appear to result from the 
adiabatic inner wall of the analytical model versus the solid 
conductor-rod of the numerical model, as the rod also showed no 
change in temperature over 40 minutes. A more direct comparison 
would result for an aspect ratio of, say, 10 in the numerical 
model, together with a program modification allowing an adiabatic 
inner wall or a thermal mass in the analytical model. 

Finally it is noted that, for comparison with the flight 
experiments, the analytical results are only as accurate as are^ 
the knowledge of the material properties. This is l 1 lustrated an 
Fig. 3-7 where the annulus surface temperature variations are 
shown parametrically with the canister surface absorptivity. It 
is seen that there can be substantial changes in the outer wall 
temperatures, resulting from changes in the canister surface 
properties. Therefore, for best comparison, it is highly 
desirable that these properties be known in the experiments. 
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<D <3 0) X 


NOMENCLATURE 


a diffusivity (cm /s) 

a,b fitted constants 2 

B ’ blackbody radiosity ( W/cm ) 2 

C modified blackbody radiosity (W/cm ) 

c radiation constant, specific heat (J/gK) 

c speed of light in vacuum (m/s) 

° emmissive power (W/cm ) ' 

emissivity fractional function 
mass velocity (g/s-cm ) 
gravitational acceleration (m/s ) 

latent heat of fusion (J/g) ru/cm 2 K) 

h Planck's constant, heat transfer coefficient ™ ) 

k extinction coefficient, thermal conductivity (U/cmK) 

N radiation Biot number 

n refractive index 

0 heat transfer rate per unit length (W/cm) 

q heat flux (W/cm 2 ) 

irradiance (U 3 /cm ) 

radiosity (W/cm 2 ) _ _ 2 

? radiation network resistance (cm ) 

r electrical resistivity (Ohm-cm), radaus (cm) 

r absolute temperature (K) 

t time (s) 

Li,v velocities (cm/s) 


q- 

q+ 


Greek 

i*y absorptivity 

' t , density (9/cm 3 ), reflectivity 

C emissivity 

r integration variable 

wave length 

0 nondimensional surface temperature 

•t nondimensional time 

f radius ratio 2 4 \ 

„ Stef an-Boltzmann constant (W/cm K ) 

r volumetric expansion coefficient (K ) 


Su bscript s 

0 pertaining to Tq 

a pertaining to transparent upper 

b blackbody 

e effective 


limit , 


1 


£ surface i 

£j surface i to surface j 

k conduction 

j> liquid 

m melting or phase-change temperature 

n normal direction 

r radiation, reference 

s solid, source 

\j V3CUUfTl 

y spectral quantity 
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Fig. 2-5. Square of Effective Refractive Index for Lithium 
Flouride versus Temperature. 
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